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Abstract 

We present a fit to precision electroweak data in the standard model extended by an additional 
vector boson, Z', with suppressed couplings to the electron compared to the Z boson, with couplings 
to the b-quark, and with mass close to the mass of the Z boson. This scenario provides an excellent 
fit to forward-backward asymmetry of the b-quark measured on the Z-pole and ±2 GeV off the 
Z-pole, and to lepton asymmetry, A e , obtained from the measurement of left-right asymmetry 
for hadronic final states, and thus it removes the tension in the determination of the weak mixing 
angle from these two measurements. It also leads to a significant improvement in the total hadronic 
cross section on the Z-pole and Rb measured at energies above the Z-pole. We explore in detail 
properties of the Z' needed to explain the data and present a model for Z' with required couplings. 
The model preserves standard model Yukawa couplings, it is anomaly free and can be embedded 
into grand unified theories. It allows a choice of parameters that does not generate any flavor 
violating couplings of the Z' to standard model fermions. Out of standard model couplings, it only 
negligibly modifies the left-handed bottom quark coupling to the Z boson and the 3rd column of 
the CKM matrix. Modifications of standard model couplings in the charged lepton sector are also 
negligible. It predicts an additional down type quark, D, with mass in a few hundred GeV range, 
and an extra lepton doublet, L, possibly much heavier than the D quark. We discuss signatures 
of the Z' at the Large Hadron Collider and calculate the Z'b production cross section which is the 
dominant production mechanism for the Z'. 
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I. INTRODUCTION 



Among the largest deviations from predictions of the standard model (SM) is the dis- 
crepancy in determination of the weak mixing angle from the LEP measurement of the 
forward-backward asymmetry of the b-quark, A b FB , and from the SLD measurement of left- 
right asymmetry for hadronic final states, y4 e (LR — had.). These two measurements, showing 
the largest deviations from SM predictions among Z-pole observables, create a very puzzling 
situation [I], [2J. Varying SM input parameters, especially the Higgs boson mass, one can 
fit the experimental value for one of them only at the expense of increasing the discrepancy 
in the other one. While A b FB prefers a heavy Higgs boson, ~ 400 GeV, A e (LR — had.) 
prefers ~ 40 GeV. Since other observables also prefer a lighter Higgs the focus has been 
on possible new physics effects that modify A b FB . However, if the pull for a large Higgs 
mass from A b FB is removed, the global fit preference is in tension with LEP exclusion limit, 
m h > 114 GeV [3]. 

In a previous study jl] we showed that a Z' with mass close to the mass of the Z boson, 
with suppressed couplings to the electron compared to the Z-boson, and with couplings to 
the b-quark, provides an excellent fit to measurements of A FB on and near the Z-pole, and 
simultaneously to A e (LR — had.). It also leads to a significant improvement in the total 
hadronic cross section on the Z-pole and Rb measured at energies above the Z-pole. In 
addition, with a proper mass, the Z' can explain the 2.3a excess of Zbb events at LEP in the 
90 — 105 GeV region of the bb invariant mass, thus expanding the family of possible expla- 
nations of the excess that include a Higgs boson with reduced coupling to the Z boson [5H7j 
or a SM-like Higgs boson with reduced branching fraction to bb [8-1 0J. 

In this paper, we explore in detail properties of the Z' needed to explain the data and 
present a model for Z' with required couplings. We discuss signatures of the model at the 
Large Hadron Collider. We calculate the Z'b production cross section which is the dominant 
production mechanism for the Z' and discuss signatures of extra vector-like quarks that are 
predicted by the model. 

We consider a new vector boson, Z', associated with a new gauge symmetry U(l)', with 
couplings to the electron and the b-quark: 

£ D Z'^(g'iP L + g'nlP R )e + Z'^{g' b P L + g' b P R )b. (1) 
In the numerical analysis we do not make any assumptions about the origin of the Z' and 
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treat all four couplings and the mass of the Z' as free parameters. Couplings to other SM 
fermions and the mixing with the Z boson are assumed to be negligible and are set to zero 
for simplicity. Once we determine typical sizes of g'^' R couplings required we construct a 
model that generates them through mixing of standard model fermions with extra vector- 
like fermions charged under the U(l)'. This method of effectively generating Z' couplings 
was recently discussed in Ref. [TTj . Although this is not the only possible model, it is a 
simple one that preserves standard model Yukawa couplings, it is anomaly free and can be 
embedded into grand unified theories. 

Previous explanations of the deviation in A b FB focused on modifying g R . Achieving this 
and simultaneously not upsetting quite precise agreement in Rf, turned out to be very chal- 
lenging for a new physics that enters through loop corrections [12]. This motivated tree level 
modification of the g R either through mixing of b-quark with extra fermions [T3J or through 
Z-Z' mixing [TH [15]. However the A b FB is only a part of the puzzle and any new physics 
that reduces to modification of bottom quark couplings cannot affect the j4 e (LR — had.). 

In Ref. [1] we suggested to modify the bb production cross section directly, e + e~ —> Z'* — > 
bb, rather than modifying the Z-couplings. This idea comes from a simple observation that 
increasing cr(ei,eL — > 6r5r) can decrease A b FB and simultaneously increase v4 e (LR — had.) by 
a relative factor that is needed to bring them close to observed values while still improving 
on Rj, (for more details see Ref. p[J). This can be only achieved by a Z' near the Z-pole. To 
generate a sizable contribution to A b FB while contributing to R^ only negligibly on the Z- 
pole, and not significantly affect predictions for A b FB and Rf, above the Z-pole (that roughly 
agree with measurements), the increase in a(eiei — > &r&r) must be due to the s-channel 
exchange of a new vector particle with mass close to the mass of the Z boson. A scalar 
particle near the Z-pole can modify A b FB only comparably to its modification of Rb- 1 A 
heavy particle, or a particle contributing in t-channel, can modify Z-pole observables only 
negligibly if it should not dramatically alter predictions above the Z-pole. Thus a Z' near 
the Z-pole with small couplings to the electron (in order to satisfy limits from searches for 
Z') and sizable couplings to the bottom quark is the only candidate. 

There is extensive literature concerning models for Z' and their phenomenological im- 

1 This was considered in Ref. |16j motivated by previous discrepancies in .Z-pole observables, namely a large 
deviation in Rb (which currently agrees with the SM prediction). 
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plications [T7j. A Z' was frequently used to explain previous discrepancies in precision 
electroweak data, see e.g. a heavy Z' [IB] or almost degenerate Z and Z' [T3] scenarios. 
Related constraints on a Z' near the Z-pole were discussed in Refs. [201 EI]- 

This paper is organized as follows. In Sec. II we outline the numerical analysis. The 
results of the best fit to precision electroweak data and ranges of Z' mass and couplings 
needed to fit the data are presented in Sec. III. A possible model leading to required couplings 
is discussed in Sec. IV. The current constraints and LHC predictions are summarized in 
Sec. V, and we conclude in Sec. VI. 

II. NUMERICAL ANALYSIS 

We construct a x 2 function of relevant quantities related to the bottom quark and elec- 
tron measured at and near the Z-pole which are summarized in Table [Tj Their precise 
definition can be found in the EWWG review (22] from which we also take the correspond- 
ing experimental values. Instead of the pole forward-backward asymmetry of the b-quark, 
A°p b B , we include three measurements of the asymmetry, at the peak and ±2 GeV from 
the peak. These are more relevant because the presence of a Z' near the Z-pole changes 
the energy dependence of the asymmetry. In addition, about 25% of the deviation in the 
pole asymmetry comes from the measurement at +2 GeV from the peak. Corresponding 
LEP averages for Rb at ±2 GeV from the peak do not exist. These are available only from 
DELPHI [23] and although they are included in the Z-pole LEP average, R®, we include 
them in addition in order to constrain the energy dependence. We further include pole 
values of the total hadronic cross section, cr° ad , the ratio of the hadronic and electron de- 
cay widths, R° e , forward-backward asymmetry of the electron, A°p B , measured at LEP; and 
SLD values of asymmetry parameters of the b-quark, A^, obtained from the measurement of 
left-right forward-backward asymmetry, and the electron, obtained from the measurement 
of left-right asymmetry for hadronic final states, A e (LR — had.), and leptonic final states, 
A e (LR — lept.). 

Although we do not modify production cross sections for the charm quark and other 
charged leptons we nevertheless include related observables in the x 2 because of correlations 
with observables related to the bottom quark and the electron. The correlations are included 
for the following two sets of observables. The first set consists of 9 pseudo-observables: 



4 



raz, Tz, o^ad? R° e i -^n Ap Bl A FB) A FB . The second set represents 18 heavy-flavor 
observables: R° b , R° c , A b FB {-2), A FB (-2), A b FB (pk), A c FB (pk), A» FB {+2), A c FB {+2), A b , A c , 
and 8 additional b- and c-tag related observables that we fix to the fitted values. Precise 
definition of these observables can be found in the EWWG review [22] from which we also 
take the corresponding experimental values and correlations. 

While b-quark quantities were measured at three energies near the Z-pole, the total 
hadronic cross section was measured also at ±1,3 GeV (from data collected only during 
1990-1991) by four LEP collaborations [24T - I27] . Because there are no combined results, we 
take ALEPH results to estimate the relative errors for each measured o"had{y/s) as 1.8%, 
0.4%, 1.1%, 1.1%, 0.3%, 1.3% for = -3, -2, -1, 1, 2, 3 GeV from the peak respectively. 2 
We then require the total hadronic cross section including Z' to deviate from the SM cross 
section no more than twice the estimated experimental error at a given yfs. 

We calculate theoretical predictions using ZFITTER 6.43 [2EJ |29] and ZEFIT 6.10 [30] 
which we modified for a Z' with free couplings to the b-quark and the electron. In the case of 
the standard model we precisely reproduce the result in the EWWG review [22j or the PDG 
review [2] for sets of SM input parameters used in those fits. In our fit we use the SM input 
parameters summarized in Table 8.1 of the EWWG review [22], namely: rriz = 91.1875 
GeV, Aa( 5 )(m|) = 0.02758, 

a s{ m z) = 0.118; however, we update the top quark mass to 
the Tevatron average, m t = 173.3 GeV [31], and fix the Higgs mass to ran = 117 GeV. As 
a result of the different set of SM input parameters, our SM predictions, given in Table [I], 
are slightly different from [21 122]- The effects of varying input parameters on electroweak 
observables can be found in [22] • The differences resulting from a given choice of SM input 
parameters are not essential for comparison of the SM and SM+Z'. We minimize the x 2 
function of 5 parameters, m z >, g'£, g% g'l, and g'^, with MINUET [52]. In principle, the width, 
V z'-, could be treated as a free parameter because Z' can have additional couplings that do 
not affect precision electroweak data. For simplicity, we do not consider this possibility. 

2 ALEPH collaboration quotes both statistical and systematic errors and the combined errors are compa- 
rable to just statistical errors quoted by other LEP collaborations. 
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III. RESULTS 



The best fits to precision data included in the x 2 are summarized in Table [I] and param- 
eters for which the best fits are obtained are given in the caption. The fit I corresponds 
to all four couplings in Eq. allowed to vary, while the fit II assumes that only two cou- 
plings, g'l and g^, are free parameters, and g% = g'l = 0. Best fits are also compared with 
predictions of the standard model. Clearly, addition of a Z' provides an excellent fit to 
selected precision electroweak data with x 2 — 6.8 with 5 additional parameters compared 
to the standard model that has x 2 = 25. 3 The most significant improvement comes from 
the three measurements of A b FB which can be fit basically at central values in the Z' model, 
without spoiling the agreement in Rb. The energy dependence of both quantities near the 
Z-pole for both the SM and the Z' model together with data points is plotted in Fig. [T} 
The A e (LR — had.) and <r° ad are also fit close to their central values. The fit II illustrates 
that most of the improvement originates from two couplings g'l and g' B which was already 
discussed in Ref. [3]. Allowing all four couplings further improves the fit and also enlarges 
the ranges of couplings for which a good fit is achieved. 

Besides quantities included in the x 2 an d given in Table [i] we check all other electroweak 
data on and near the Z-pole, and above and below the Z-pole. While b-quark quantities were 
measured at three energies near the Z-pole, the total hadronic cross section was measured 
also at ±1,3 GeV as discussed in the previous section. The measurement at +1 GeV 
roughly coincides with the Z'-peak where the deviation from the SM would be the largest. 
The experimental error in cihad at +1 GeV from the peak is ~ 1% for each LEP experiment 
and thus the Z'-peak contributes only a fraction of the error bar. 

At energies above the Z-pole, the A b FB in the Z' model basically coincides with the SM 
prediction while Rb fits data better than the SM, see Fig. [TJ with x 2 — 4.9 for 10 data points 
compared to the SM which has x 2 = 7.3 (the average discrepancy with respect to the SM 
prediction for R b is —2.1a) |33j. At energies below the Z-pole the Z' leads only to negligible 
differences from the SM predictions compared to sensitivities of current experiments. 

The quantities related to other charged leptons and quarks are not directly affected by Z' 

3 The difference in % 2 compared to Ref. [3] is the result of a more complete \ 2 function that includes more 
observables and correlations between them. The best fit values of Z' parameters and main features of 
results are not affected by these modifications. 
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TABLE I: The best fits to relevant precision electroweak observables in the SM with a Z' . The 
best fit I assumes all 4 couplings are free parameters and it is achieved for: mz' = 92.2 GeV, 
g'l = 0.0065, g% = 0.0077, g'£ = 0.044, and g% = -0.51; (I> = 1.0 GeV). The best fit II 
assumes only two couplings are free parameters (g'^ = g'l = 0) and it is achieved for: mz> = 92.2 
GeV, g'l = 0.0048, and g'^ = —0.52; {T z 1 = 1-0 GeV). The standard model input parameters 
are fixed to rrit = 173.3 GeV, rrih = 117 GeV, and other parameters as listed in Table 8.1 of the 
EWWG review [22] . For comparison, we also include predictions of the standard model with x 2 
contributions. 



Quantity 


Exp. value 


SM 


2 

XSM 


I 


2 

Xi 


II 


Xii 


D r i l 
^had i nb ] 


41.541(37) 


41.481 


2.6 


41.541 


0.0 


41.532 


0.1 


R b (-2) 


0.2142(27) 


0.2150 


0.1 


0.2157 


0.3 


0.2161 


0.5 


R° b 


0.21629(66) 


0.21580 


0.6 


0.21693 


1.0 


0.21676 


0.5 


R b {+2) 


0.2177(24) 


0.2155 


0.8 


0.2179 


0.0 


0.2185 


0.1 


A b FB (-2) 


0.0560(66) 


0.0638 


1.4 


0.0581 


0.1 


0.0583 


0.1 


A b FB (pk) 


0.0982(17) 


0.1014 


3.5 


0.0980 


0.0 


0.0971 


0.4 


A h FB {+2) 


0.1125(55) 


0.1255 


5.6 


0.1133 


0.0 


0.1139 


0.1 


A h 


0.924(20) 


0.935 


0.3 


0.921 


0.0 


0.923 


0.0 


R° e 


20.804(50) 


20.737 


1.8 


20.772 


0.4 


20.759 


0.8 


,0,e 
FB 


0.0145(25) 


0.0165 


0.7 


0.0176 


1.6 


0.0167 


0.7 


,4 e (LR-had.) 


0.15138(216) 


0.14739 


3.4 


0.15047 


0.2 


0.14849 


1.8 


A e (LR-lept.) 


0.1544(60) 


0.1473 


1.4 


0.1473 


1.4 


0.1474 


1.4 


total x 2 






24.6 




6.76 




9.99 



and the predictions are essentially identical to predictions of the SM [2]. For example, the 
LEP 1 average of leptonic asymmetry assuming lepton universality, A\ = 0.1481 ± 0.0027, 
agrees very well with the SM prediction and would be only negligibly altered by the Z' with 
couplings corresponding to the best fit (the prediction is the same as for v4 e (LR — lept.) 
given m Table (T)). 

The x 2 is n °t very sensitive to exact values of couplings. As can be seen from contours 
of constant x 2 m 9l ~ 9r P iane given in Fig. [2J a significantly better fit compared to the 
standard model can be achieved in a large range of couplings. Contours of the other two 
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FIG. 1: Experimental values of A b FB (top) and Rb (bottom) and predictions of the SM (thin lines) 
and the Z' model (thick lines) for input parameters corresponding to the best fit I specified in the 
caption of Table [I] as functions of center of mass energy near and above the Z-pole. 



couplings, g'x and g F) corresponding to the best fit in the g'£ — g% plane are given in Fig. |3 
and contours of constant mz> and the width of Z 1 determined from its couplings to the 
electron and the bottom quark are given in Fig. |4j In these and following plots the x 2 
contours from Fig. [2] are overlayed for easy guidance of the fit quality. 

Partial contributions to \ 2 from each observable are given in Figs. [5] - m Those from 
Table [i] that are not plotted, namely A&, A° F e B , and A e (LR — lept.), vary negligibly with 
varying the couplings. From these plots we clearly see that the main drivers toward the 
region of the best fit are A b FB (+2), given in Fig. [6] (Right), disfavoring the upper- left region 
of couplings in the g'£ — g'^ plane, and A e (LR — had.), given in Fig. [7] (Middle), disfavoring 
the lower-right region. The three measurements of Rb further constrain the allowed region 
of couplings approximately along the diagonal, see Fig. [5} The cx° ad , given in Fig. [7] (Left), 
fits close to the central value in a large range of couplings, and finally, R° e prefers the central 
and lower region in the g'£ — g' B plane, see Fig. [7] (Right). 

We have seen in Fig. H (Left) that the minimum of x 2 prefers close to the mass of the 
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FIG. 2: Contours of constant x 2 • Couplings g'£ and g'£ are free parameters. The shaded region 
corresponds to x 2 > 15. 
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FIG. 3: Contours of constant g B (Left) and g'£ (Right) from the best fit in the g'£ — g B plane with 
X 2 contours from Fig. [2] overlayed. 

Z boson with the best fit requiring Z' just ~ 1 GeV heavier than the Z boson, mz> = 92.2 
GeV. However, much better fit compared to the standard model can be obtained even for 
somewhat heavier Z'. Minimum of x 2 as a function of mz> for the fit with all for couplings 
being free parameters and for the fit with only two couplings being free parameters is plotted 
in Fig. [8j In the same figure, we also show the best fit with A b FB (+2) removed from the x 2 
function near the region of the best fit. This x 2 function is almost flat which demonstrates 
that the best fit value of mz> is mainly driven by the +2 GeV measurement of the A FB . 

Finally, let us comment on the case with only 2 allowed couplings when moving away 
from the best fit presented in Table [ij Contours of constant x 2 m 9l ~ 9r plaue in the case 
when only these two couplings are allowed, and thus <7jf = g% — 0, are given in Fig. [9] (Left). 
Preferred region of these couplings is very similar to the case with all four couplings allowed, 
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FIG. 4: Contours of constant mz' (GeV) (Left) and the width of Z 1 (GeV) (Right) with x 2 contours 
from Fig. [2] overlayed. 
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FIG. 5: Contours of constant contribution to x 2 from i?&(— 2) (Left), (Middle), and R\,{+2) 
(Right) with x 2 contours from Fig. [2] overlayed. 
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FIG. 6: Contours of constant contribution to x 2 horn A FB (—2) (Left), A b FB (pk) (Middle), and 
A b FB (+2) (Right) with x 2 contours from Fig. [2] overlayed. 
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FIG. 7: Contours of constant contribution to x 2 from cr° ad (Left), A e (LR — had.) (Middle), and 
i?g (Right) with x 2 contours from Fig. [2] overlayed. 
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FIG. 8: Minimum of x 2 as a function of mz' for the fit with all for couplings being free parameters 
(middle line), and for the fit with only two couplings being free parameters, assuming g B = g F = 
(top line). In addition, a fit without A b FB (+2) in the x 2 function near the region of the best fit is 
shown (bottom line) demonstrating that the best fit value of mz> is mainly driven by the +2 GeV 
measurement of the A h FB . 



see Fig. [2J which demonstrates that g'£ and g'^ are the relevant couplings responsible for 
dramatic improvement of the fit compared to the standard model. Contours of constant 
mz'y given in Fig. [9] (Right), also closely resemble those of the four coupling fit, see Fig. [4 
(Left). The main difference from the previous fit with all four couplings allowed and the 



reason for somewhat worse \ 2 are ^fbIp^) and A e {LK— had.) given in Fig. 10 Contributions 
to x 2 from other observables is very similar to the previous case with all four couplings. 
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FIG. 9: Contours of constant x 2 (Left) and mz> (GeV) (solid lines) with x 2 contours (dashed lines) 
(Right) with g'ft = g'£ = 0. The shaded region corresponds to x 2 > 15. 
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FIG. 10: Contours of constant contribution to x 2 from ^4^ s (pk) (Left) and A e (LH — had.) (Right) 
with g'ft = g'£ = 0. The x 2 contours from Fig. [9] are overlayed. 

IV. A POSSIBLE MODEL AND ITS CONSEQUENCES 

In order to have large enough contribution to Z-pole observables without significantly 
modifying above the Z-pole measurements the mass of the Z' should be within few GeV 
from the Z mass. Couplings that are required are: g'[ ~ 0.005 and g% ~ —0.5. Additional 
small g'£ and g'l further improve the fit to Z-pole data but are not required. Their presence 
however expands regions of g'£ and g'^ where a good fit is achieved. 

The required couplings of standard model fermions to Z' do not follow the usual pattern 
expected from new gauge interactions. A simple framework to generate arbitrary couplings 
of standard model fermions to a Z' while preserving Yukawa interactions and keeping the 
model anomaly free was recently discussed in Ref. [11]. In this framework the couplings of 
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standard model fields to Z' are generated effectively through mixing with extra vector-like 
fermion pairs. We will follow this direction and customize it for our purposes. 



TABLE II: Quantum numbers of relevant standard model and extra vector-like particles. 

QL d R l L e R H D L D R L l L r $ 

SU(3) C 33111 3 31 11 

SU(2) L 21212 11221 
tt/i \ l l l 1 l l l l l n 

U(l)' -I -I 1 I -1 



Let us start with adding a vector-like pair of fermions Dl and D R charged under a U(l)' 
where D R has the same quantum numbers under the standard model gauge symmetry as 
the d R , see Table [XTJ This charge assignment results in the following renormalizable terms 
in the lagrangian: 

C D -q u Y£d Rj H - D L \ d k d Rk $ - fi D D L D R + h.c. (2) 

where the first term represents the usual standard model Yukawa couplings for down-type 
quarks (the sum over flavor indices is assumed). The second term contains Yukawa inter- 
actions of standard model quarks and the extra Dl quark. The last term is the mass term 
for the vector-like pair. The vacuum expectation value of the scalar field $ breaks the U(l)' 
and generates mixing terms between d R i and D R . After spontaneous symmetry breaking 
the 4x4 mass matrix for down type quarks is given by: 

\D R ) \\j(<I>) no) \D R ) 

and it can be diagonalized by a bi-unitary transformation, U L MdU R , which defines the mass 
eigenstate basis. However, before doing that, it is instructive to change the basis by a unitary 
transformation, d R —> V R d R , d L — > Vid L) which diagonalizes the standard model Yukawa 
couplings, Y d . The mass matrix becomes: 

\ d n V Rnj (<f>) n D ) \aj l) 
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where 

^w*>_ wd (vjy^m. (5) 

From this form of the mass matrix we can see that in any theory of flavor that determines the 
structure of Yukawa matrices for standard model fermions, in this case only Yd, but allows 
arbitrary Xf couplings, these can be chosen so that a± = «2 = and only 0:3 is non-zero. 
This corresponds to the situation when Xf oc V R3i , or in the basis where standard model 
Yukawa couplings are diagonal, it corresponds to Xf = Xf = and A3 = A& is non-zero. This 
is the minimal scenario that does not modify standard model couplings of down and strange 
quarks. In what follows we will focus on this scenario. 

In the basis where standard model Yukawa couplings are diagonal, assuming Xf are such 
that q;i = «2 = 0, the first two diagonal entries correspond to masses of the down and 
strange quarks: 

m d>s = n D /3i ;2 (6) 

The lower 2x2 block can be diagonalized by a bi-unitary transformation (for simplicity we 
drop indices, a 3 = a and (3 3 = (3): 



HoUl r '\U R =\ " I , (7) 



where we use the same names, Ul,r, for matrices that diagonalize the lower 2x2 block in the 
case ot\ t 2 = 0, as for the matrices that diagonalize the general 4x4 matrix. We label their 
components by 3 and 4 so that results are applicable to the general scenario with non-zero 
ax 2- The bottom quark mass and the mass of the extra heavy down-type quark are given 
by: 



1/3 0\ 




' m h \ 




[a 11 




1 m D j 



m b ~ fi D p/Vl + a 2 , (8) 
m D ~ hdVI + a 2 , (9) 

where we assume /3<l,a. The mass of the D quark as a function of /xp and a is plotted 



in Fig. 11 The diagonalization matrices are approximately given by: 




"1* u R ~-^=\ - |. (10) 
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500 




FIG. 11: The mass of the extra D quark, mo [GeV], as a function of hd and a. 
A. Couplings of the Z' boson 

Couplings of Z' to down-type quarks (mass eigenstates) originate from the kinetic term 
of the extra vector-like pair: 

C ki n D D L 0'D L + DrHJ)'D R = Li{Ul)iA0'{U L ) Aj d Lj + Li(U R )i4i$'(U R ) 4j d Rj , (11) 

where the vectors of mass eigenstates are da = (d R , s,r,&r, Dr) t and similarly for d^. The 
covariant derivative is given by: 



D'» = D™ - ig'Q'Z'^ 

where Df, M is the standard model covariant derivative: 

9 



D SM = a _■ 

» COS0 



W 



(T 3 -sin 2 6 w Q)Z IM -ieQA li , 



(12) 



(13) 



and for simplicity we do not write the SU(3)c interactions explicitly which are not modified 
by field redefinitions. In the mass eigenstate basis, the Z' has in general both flavor diagonal 
and off-diagonal couplings to down-type quarks: 

9r U = -9' (U R MU R ) 4j (14) 
9l U = -9 (Ul)u(U L ) 4j: (15) 

where we used Q' D = — 1. For flavor diagonal couplings the expressions simplify to: 

9t = -9' WrU 2 (16) 
9't = -9 \{UlU 2 - (17) 
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FIG. 12: Contours of constant -g'j^ (Left) and -g'l (Right) in hd - a plane for g 1 = 1. 




FIG. 13: Contours of constant -g'j^ (Left) and -g'£ (Right) in g' - a plane for /j,d = 200 GeV. 

In the case atip = 0, that we are focusing on, the first two generations do not have couplings 
to Z' and only the bottom quark and the D quark couple to Z' with couplings that can be 



obtained from Eq. (10). The g'^ L couplings as functions of [in and a assuming g' = 1 are 



given in Fig. 12, and as functions of g' and a, for fixed = 200 GeV, in Fig. 13 The 
g% coupling is fully controlled by a and can be easily sizable. For g' = 1 it can be close to 
the value suggested by the best fit (highlighted in plots) for a ~ 1. The g'l coupling on the 
other hand is proportional to (3 which is of order m b / jj D , see Eq. (pi), and thus it is very 



small. For the purposes of the fit to precision electroweak data it is effectively zero. 
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B. Corrections to neutral and charged currents 



All couplings of the photon and couplings of up-type quarks and right-handed down-type 
quarks to the Z boson are identical to standard model couplings. However, since Dl is an 
SU(2)l singlet, the couplings of left-handed down-type quarks to the Z boson are modified. 



They can be read out from kinetic terms of left handed fields (similar to Eq. ( 11 ) but written 
for all four quarks): 

4 

& U = -^Y,^-^ 2e wQd)iUl) tk {U L ) k ^ (18) 
cos U w ^ 

where T| = —1/2 for k = 1, 2, 3 and for k — 4. Corrections to couplings of the Z boson 
to left-handed down-type quarks in the standard model (i, j = 1, 2, 3) can be written as: 

Sg f L ifi = i^r( u lU u L)^ (19) 

Z COS U\y 

In general, these corrections for the first two generations are tiny, since they are proportional 

f f ■ 

to ratios of masses of corresponding quarks and the heavy quark, 5g£ 3 oc (m^/ nr,){mj/ ^d)- 
In the case a± t 2 = 0, that we are focusing on, couplings of the first two generations to Z are 
not altered at all, and there are no flavor violating couplings. 



Comparing Eq. (19) with Eqs. (15) and (17) we see that the change in a Z coupling is 



directly proportional to corresponding Z' coupling that is being generated. For the correction 
to the left-handed bottom coupling we find: 

Sg b L = -^r4- (2°) 

2costV g 



and from the values of the ratio g'l/g' given in Fig. 12 we see that Sg b L is negligible. 
The charge currents, 

- -^u Li {VcKM)irfd Lj W+ + h.c. (21) 
get also modified by — > U id^ which effectively leads to a modification of the CKM matrix: 

3 

{VcKM)ij -> ^ V CKM)ik{U L ) kj . (22) 
k=l 

In the case a± 2 = 0, only the third column of the matrix is modified: 

(VcKM)ib (VcKAl)ib(U L )33, i = u,c,t. (23) 
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FIG. 14: Relative correction 5 to the third column of the CKM matrix. 
It is convenient to define 5 = 1 — (Ul)^ which represents the relative correction of the third 



column of the CKM matrix. It is plotted in Fig. [14] and the values are far below current 
uncertainties in the CKM elements. 



C. Exploring lagrangian parameters 

The model we have discussed so far is specified by 4 parameters: g', A&, //£>, and the 
vacuum expectation value (vev) of the extra Higgs field that breaks the U(l)' symmetry, 
($). This vev is responsible for generating the Z' couplings to b quark through mixing with 
D (it is contained in a) and also for the mass term of the Z' boson: 

D \D^\ 2 D g' 2 {^) 2 Z[ L Z^, (24) 



m z > = y/2g'($). (25) 

Equivalently, the model is specified by g', mz', a, and although the fit to precision 
electroweak data depends only on two parameters: mz> and g'^. The fit strongly prefers 
mass of the Z' close to the mass of the Z boson and thus we can simply fix it to the best 
fit value 92.2 GeV. In previous subsections we have explored the dependence of g^ coupling 
on g', a, and /ip. It is however instructive to see what values of lagrangian parameters and 
($) are required. Contours of constant g% in the g' - Xb plane for values of \id = 100, 200, 



and 500 GeV are given in Fig. 15 The corresponding vev of $ is given on the right axis 
and the mass of the D quark is overlayed. In order to obtain g^ ~ 0.5, as suggested by the 
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FIG. 15: Contours of constant g'^ (solid lines) in the g' - Xb plane for fi£, = 100 GeV (Top Left), 
200 GeV (Top Right), and 500 GeV (Bottom). The mass of the Z' is fixed to the best fit value, 
mz' = 92.2 GeV, and the corresponding vev of the extra Higgs field is given on the right axis. The 
dashed lines represent contours of constant D quark mass, mo [GeV] (vertical line = which is 
not shown would correspond to mp = fio)- 

best fit, while keeping g' and Xb perturbative, the extra D quark should be fairly light, in 
a few hundred GeV range. However we should keep in mind that even ~ 0.1 provides 
a significant improvement of the fit compared to the standard model, in which case the D 
quark can be heavier. 

For completeness, we also plot contours of constant g'^ D in the g' - Xb plane for /xp = 200 



GeV in Fig. 16 
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FIG. 16: Contours of g'^ D (solid lines) in the g' - Aft plane for /j,£> = 200 GeV. The rest as in Fig. 15 
D. Coupling of the electron to Z' 

The other required coupling besides g'^ is g'£. This coupling can be generated in a very 
similar way, by adding a vector-like pair of fermions Ll and Lr charged under the U(l)' 
where Ll has the same quantum numbers under the standard model gauge symmetry as 
the lepton doublet II, see Table [XTJ The U(l)' charge assignment for heavy fermions and 
$ is chosen so that heavy fermions fit into complete GUT multiplets, in this case 5 and 5 
of 577(5). The model is anomaly free and its supersymmetric version preserves the gauge 
coupling unification of the standard model gauge couplings. 

Our charge assignment results in the following renormalizable terms in the lagrangian: 



C D -luY^eRjH - l Lk \ l k L R & - HlL l L r + h.c. 



(26) 



where the first term represents the usual standard model Yukawa couplings for charged 
leptons. The second term contains Yukawa interaction between lepton doublets and the 
extra Lr lepton. The last term is the mass term for the vector-like pair. 

The derivation of couplings in the charged lepton sector and the discussion of flavor 
violation closely follow the down quark sector. In the basis where standard model Yukawa 
couplings are diagonal we choose A' = (A e ,0,0). This is the minimal case which generates 
g'l while the couplings of \x and r to Z' and flavor violating couplings are not generated. 

Due to opposite Q' charges of Lrx and Dr^l, motivated by an SU(5) embedding, cou- 
plings g'l and g f R have automatically opposite signs which is required by the best fit. The 
dependence of g'l on parameters of the model is identical to what we presented for g'^, how- 
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ever the value of interest is much smaller. The value of g'£ motivated by the best fit is about 
1% of the g'ft, see Table [i] This can be achieved for: 

^ ~ 10 ^, (27) 

which means that either the mixing coupling A e is very small compared to \ or the extra 
lepton L is much heavier than the extra down quark D. Since the mass of the electron is 
negligible, the generated g'£ and corrections to Z couplings are essentially zero. 



E. Extensions of the model for other Z' couplings and Z-Z' mixing 

So far we have considered a model that adds vector-like fields, with charges consistent 
with embedding into 5 and 5 of SU(5). Such a model generates only g'^ and g'£ couplings. 
The best fit with just these couplings is the fit II in Table [TJ Additional small couplings, g'£ 
and g'l, improve the quality of the fit somewhat. Generating even sizable g'£ presents no 
challenge. However, g'l leads to a modification of both the 3rd row and 3rd column of the 
CKM matrix. More importantly these corrections are not suppressed by the mass of the b 
quark and thus the generated g'l cannot be very large. However, the value of g'l suggested 
by the best fit to precision electroweak data is quite small, see Fig. [3] (Right), and even 
g'l = does not significantly change the fit. The g'l is the least important coupling of the 
four. One can consider generating these additional small couplings by adding vector-like 
fields with charges consistent with embedding into 10 and 10 of SU(5). 

Since the extra vector-like fermions couple to both Z and Z' their loops can generate 
Z-Z' mixing. The contribution of a vector-like pair to the mixing can be however cancelled 
by adding a second vector-like pair with opposite U(l)' charge. In addition, the mixing can 
be avoided when the U(l)' is embedded into a non-abelian group. 



V. Z' AND D AT THE LARGE HADRON COLLIDER 



At hadron colliders the Z' could be produced in association with b quarks, See Fig. 17 



The production cross sections of Z'b at the LHC are shown in Fig. 18 for center-of-mass 
energy of 7 TeV (Left) and 14 TeV (Right). The cross sections are calculated with MCFM 
at the leading order (LO). We used CTEQ6.6 parton distribution functions (PDF). The 
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FIG. 17: Feynman diagram for Z' production in association with the b quark. 
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FIG. 18: Z'b production cross section (nb) at the LHC with center-of-mass energy 7 TeV (Left) 
and 14 TeV (Right) with x 2 contours from Fig.[2]overlayed. The shaded upper regions are excluded 
by the CDF search for the Higgs boson assuming B(Z' — > bb) = 100%. 

factorization and renormalization scales are set to fip = = Mz- For the final state b-jet, 
p h T > 15 GeV, |7;| < 2.5 and AR < 0.7 are chosen to match those used for the calculation 
of Zb production which is a background for Higgs searches [35] . In the analysis no tagging 
efficiencies are assumed. 

From Fig. 18 we see that the cross section is only governed by g'^ since g'l is negligibly 



small. In the region of the best fit the Z'b cross sections are ~0.5 nb at 7 TeV and ~2 nb at 
14 TeV. If other couplings besides g'[' R are absent the Z' would decay to bb with branching 
ratio close to 100%. The search for the Z' is therefore very similar to the search for the 
Higgs boson in 3b final state. 

Recent limits on a(pp — > <pb) x BR((p — » bb) set by CDF with 2.6 fb _1 of integrated 
luminosity [35] and DO with 5.2 fb _1 [37J already constrain the allowed values of g^ . We 
calculated the production cross sections of Z'b at the LO using MCFM with the center-of- 
mass energy of 1.96 TeV, p h T > 15 GeV, \rj\ < 2 and AR < 0.4 that are used in the CDF 
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search which currently gives strongest limits. Comparing it with the CDF limit a(pp — > 
(f)b) x BR(<f) — > bb) < 26.4 pb for = 90 GeV, we find that g'J^ larger than 0.56 is excluded 



as shown in Fig. 18 Note however that with possible couplings of Z' to other quarks (or 
particles beyond the SM) the BR(Z' — > bb) can be highly reduced resulting in weaker limits. 

At the LHC the Z'b cross section is two orders of magnitude larger than at the Tevatron. 
So it is just a question of accumulating enough luminosity to see the signal of Z' . A search 
for the Higgs boson produced in association with the b quark has not been performed yet 
at the LHC. Since predictions for the production cross sections depend on cuts used in an 
analysis, let us make few comments. In the recent ATLAS measurement of the cross section 
for b-jets produced in association with a Z boson decaying into two charged leptons, a b-jet 
was identified with p h T > 25 GeV and \y\ < 2.1 [55] . With these cuts on px and the Z'b 



production cross section is reduced to about half of those given in Fig. [18} Note also that 
MCFM is not interfaced to parton shower/hadronisation fragmentation package, and it does 
not include multiple parton interaction (MPI). We expect about 10% change in the cross 



sections given in Fig. 18 once those corrections are taken into account |3B]- At the same 
time, the uncertainties stemming from the next-to-LO calculation, the scale dependences, 
PDF, and a s are expected to be 20%, 10%, 3%, and 2%, respectively [38] . 

The model discussed in the previous section predicts the extra D quark in a few hundred 
GeV range. Constraints from searches for the 4th generation do not apply, since the D 
quark decays into Z'b with Z' — > bb. At the LHC the D quark can be pair produced by 
QCD interactions leading to 6b final states. Since Z' — > e + e~ is suppressed compared to 
Z' — > bb by {g'[/ g'tt) 2 — 10~ 4 the e + e~4b final states are very rare. 



VI. CONCLUSIONS AND OUTLOOK 



The Z' near the Z-pole with couplings to the electron and the b-quark can resolve the 
puzzle in precision electroweak data by explaining the two largest deviations from SM pre- 
dictions among Z-pole observables: A b FB and A e (LR — had.). It nicely fits the energy de- 
pendence of A b FB near the Z-pole and improves on a{J ad on the Z-pole and Rb measured at 
energies above the Z-pole. 

We constructed a model that generates the minimal set of required couplings through 
mixing of standard model fermions with extra vector-like fermions charged under the U(l)'. 



23 



It preserves standard model Yukawa couplings, it is anomaly free and can be embedded into 
grand unified theories. The model allows a choice of parameters that does not generate any 
flavor violating couplings of the Z' to standard model fermions. Out of standard model 
couplings, it only negligibly modifies the left-handed bottom quark coupling to the Z boson 
and the 3rd column of the CKM matrix. Modifications of standard model couplings in the 
charged lepton sector are also negligible. It predicts an additional down type quark, D, with 
mass in a few hundred GeV range, and an extra lepton doublet, L, possibly much heavier 
than the D quark. 

At the LHC the Z' could be produced in association with b-quarks. The production cross 
sections of Z'b are large, in the region of the best fit as large as ~0.5 nb for center-of-mass 
energy of 7 TeV and ~2 nb for 14 TeV. If other couplings besides g'l ,b R are absent the Z' 
would decay to bb with branching ratio close to 100%. The search for the Z' is therefore very 
similar to the search for the Higgs boson in 3b final state. However, with possible couplings 
of Z' to other quarks (or particles beyond the SM) the BR(Z' — > bb) can be highly reduced 
which could make the search for Z' difficult. The optimal experiment to confirm or rule 
out the possibility of a Z' near the Z-pole would be the future linear collider, especially the 
GigaZ option, which would allow more accurate exploration of the Z-peak. 

The extra D quark can be pair produced at the LHC by QCD interactions. It dominantly 
decays into Z'b leading to 6b final states. The e + e~4b final states are highly suppressed. 

Considering other flavor conserving couplings, or small flavor violating couplings, expands 
the range of observables to which this Z' could contribute. It would be interesting to see 
if it can simultaneously explain some other deviations from SM predictions. For example 
with additional couplings in the charged lepton sector the deviation in the muon g — 2 can 
be explained [39] . However, adding any additional couplings leads to many new constraints 
that have to be carefully examined. 
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